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Chapter  I 
Introduction 

Problems  frequently  encountered  in  complex  multiple  input-ouput  systems 
stem  from  interaction  or  coupling  between  the  controlled  variables.    For  ex- 
ample, in  the  lateral  control  of  an  aircraft,  a  change  in  the  angle  of  bank  by 
aileron  deflection  will  result  in  a  change  in  yaw  angle  provided  that  there  is 
no  rudder  input  to  offset  the  change.    To  eliminate  this  interaction  between 
the  controlled  outputs,   the  theory  of  non- interacting  controls  was  developed. 
For  a  control  system  designed  in  accordance  with  this  theory,   a  desired 
change  in  one  output  variable  will  not  cause  a  change  in  the  other  output 
variables. 

The  general  method  for  designing  such  non- inter  acting  control  systems 
was  first  developed  by  A.   S.  Boksenbom  and  Ro  Hood,  Ref.  1,   and  applied 
by  them  to  the  control  of  a  turbo-jet  engine.     This  work  was  later  repeated  by 
H„  S.  Tsien,  in  ReL  2.    In  ReL  3,  H,  Freeman  extended  the  theory  to  obtain 
a  desired  degree  of  interaction  in  an  m  input  and  n  output  system.    In  Ref.  4, 
Freeman  discussed  stability  and  realizability  but  he  did  not  discuss  the  practical 
design  problems  or  use  the  more  realistic  open  loop  approach. 

This  study  will  develop  the  basic  criteria  for  non-interaction  of  the 
variables  in  an  n  input  and  n  output  system.    In  addition,   the  effect  of  dis- 
turbance inputs  and  errors  in  mathematically  describing  the  physical  system 
and  synthesizing  the  controller  will  be  discussed.    It  will  also  introduce  the 
potential  designer  of  a  non- interacting  control  system  to  some  of  the  practical 
problems  which  are  often  encountered  in  applying  the  basic  non-interaction 
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criteria.    As  a  final  step,  the  theory  will  be  applied  to  the  design  of  an  auto- 
pilot system  in  order  to  illustrate  the  design  techniques  and  the  difficulties 
that  may  be  encountered  in  a  practical  application. 


Chapter  II 
General  Theory 

In  this  chapter,  the  basic  relations  between  controller  and  controlled 
elements  of  a  system  to  effect  non-interaction  of  variables  will  be  developed. 
To  achieve  this,   a  specific  system  configuration  will  be  assumed.    However, 
the  extension  «f  this  theory  to  other  configurations  is  feasible  and  should  be 
apparent  from  the  subsequent  development. 
Introduction 

A  linear  physical  system  with  n  inputs,  d^  ,  and  n  outputs,  U-  ,  can  be 
simply  described  in  vector  notation,    In  such  notation,    J~   and   U    are  con- 
sidered as  n  dimensional  vectors  with  components  £■    and  U   .    The  coefficients 
of  the  equations  relating  w    to   &  in  the  physical  system  form  an  n  by  n 
matrix,   A    ,  which  is  composed  of  elements  a...    Note  that  the  a.,  are  oper- 
ators expressable  as  a..(p)  where  p  =  d/dt    Thus, 
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r     I  =   I  f o  n 


(2-1) 


J  =  ' 
or  in  vector-matrix  notation, 


y  =  ^  %  (2-2) 


In  a  block  diagram,   the  system  described  above  can  be  represented  as: 

I — — 


—  1 


In  the  case  where  there  are  m  inputs  and  n  outputs  and  m  is  greater  than  n, 
the  m-n  additional  inputs  are  fed  directly  through  the  system  to  become  system 


outputs.    This  implies  that  the  A  matrix  is  now  an  m  x  m  matrix  with  the 
added  CL^  elements  having  a  value  of  unity.    Such  a  case  is  described  at 
length  in  Refs.  1  and  2. 

To  achieve  non-interaction,  it  is  desired  that  u  =  x  where  X  is  an  n 
dimensional  input  vector  with  components  x..    In  designing  a  linear  system 
to  satisfy  the  above  input-output  requirements,  the  foremost  problem  is  to 
use  the  error  vector,  6=  U-X,  to  determine  ^   for  non-interaction  and  at 
the  same  time  ensure  that  6-^0.    One  possible  system  configuration  for 
achieving  these  ends  is  shown  in  Fig.   1(a). 

In  order  to  keep  the  treatment  of  the  subject  sufficiently  general,  servos, 
transducers  and  controllers  have  been  included  in  the  overall  system  shown 
in  Fig.   1(a).    The  transducers  have  been  supposed  in  light  of  the  fact  that  the 
real  outputs,   UL    ,  of  a  physical  system  can  not  be  obtained  directly,.    Thus 
the  measured  outputs,  U^    ,  are  related  to  the  real  outputs,   U[   ,  through  the 
transducer  matrix,   T.    This  matrix  is  a  diagonal  matrix  with  elements  t^. 
It  is  further  supposed  that  the  transducers  are  selected  such  that  the  Lji  —<-/<; 
are  very  small  and  hence  no  attempt  will  be  made  to  correct   U^   to  U^    . 
The  physical  system  inputs,   hi   ,  are  generally  servo  driven  and  the     d /    are 
related  to  the  servo  signal  inputs,   o'^  ,   through  the  servo  matrix,  S.     Like 
the  T  matrix,   S  is  a  diagonal  matrix  with  elements  s..o    The  controller,  which 
relates  the  system  input  signal  vector,  0    »  to  the  error  signal  vector,  £   , 
is  represented  by  an  n  by  n  matrix,   C.    It  will  be  designed  such  that  the 
overall  system  is  non- interacting  and  has  desired  response  characteristics. 

As  a  means  of  simplifying  the  ensuing  mathematical  analysis,  the  system 
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Typical  Control  System  Configurations 


shown  in  Fig.   1(b)  shall  be  used.    It  should  be  noted  that  this  is  an  equivalent 
system  to  that  of  Fig.  1(a)  and  that  many  multiple  input-output  systems  can  be 
put  into  this  form.     This  does  not  mean  that  this  is  the  only  form  possible. 
However,  conditions  for  non-interaction  are  more  complex  and  system  design 
becomes  exceedingly  difficult,  if  not  impossible,  when  configurations  different 
from  Figs.  l(af  and  1(b)  are  used.    A  simple  exception  to  the  above  is  shown  in 
Fig.  1(c).    Here  the  diagonal  matrix  R,  which  might  for  example  provide  deriva- 
tive feedbacks,  has  been  added,.    Although  another  loop  has  been  added,  the 
system  of  Fig.  1(c)  reduces  to  that  of  Fig.   1(b)  if  A  =  (T  +  R)AS  and  u  =  (T+  (?)£7. 
Analysis  of  the  Control  System 

The  system  under  consideration  is  the  simplified  but  equivalent  system 
shown  in  Fig.   1(b).    The  vector  equation  for  this  system  is: 

cj=  AS  =  ACe  =  AC(x-lj)  (2-3) 

Solving  equation  (2-3)  for  the  system  output  vector,  U   ,  provides: 

<j-     [X-*-AC]"'aC  x  (2-4) 

where  I  is  the  identity  matrix. 

In  order  that  the  inverse  matrix    [1  +  ACJ     shall  exist,  it  is  required 
that  the  determinant,  [I-I-ACJ  ,  is  not  identically  zero  over  all  values  of  the 
complex  variable  p.    Of  course,  AC  will  be  a  function  of  the  operator  p  and 
hence  jX-t-AC|  may  be  equal  to  zero  for  a  few  particular  values  of  the  operator 
p  considered  as  a  complex  variable. 

If  the  system  under  consideration  is  to  be  non- interacting  it  is  required  that: 


y 


^    DX 


(2-5) 


where  D   is  a  diagonal  matrix  with  elements  cL..    This  also  implies  that 

Ui  -  diL^-i      *or  ^  i  or  tnat  a  specific  output  is  affected  only  by  a  corresponding 

input.  * 

From  equations  (2-4)  and  (2-5)  it  follows  that: 


[ioacJ'ac  -  D 


(2-6) 


Certainly,  one  condition  on  the  matrix  AC  that  satisfies  the  above  equation  is 
that  AC  is  a  diagonal  matriXo    If  AC  is  diagonal  thenuZ-f-ACJ  AC   is  given  by: 


[i+acT'ac  = 


(ac)„ 


o 


I+(OC)m 

o    i^fo 

0 


0         0 


where  (ac)..  are  the  components  of  AC. 

Solving  equation  (2-6)  for  AC  yields: 


AC  -  [x-35]"  D 


•   •    o 
•  .  o 

o 
°  i -Hack 


(2-7) 


(2-8) 


Again,  it  is  obvious  that  AC  is  a  diagonal  matrix  if  equation  (2-8)  is  to  be 
satisfied. 

In  summarizing,   a  necessary  and  sufficient  condition  for  non-interaction 
is  that  the  matrix  AC  be  diagonal;  I.e.  AC  -  D.     Furthermore,  the  elements 


(ac)     or  d     must  be  finite  and  not  identically  minus  one.    This  assures  that: 
ii         ii 

L\l=,   dLiXi    =    ^V   Xt   «    ^liL       £rl*li»vi  (2-9) 

where  d^  is  finite. 

It  is  now  known  that  AC  =  D  where  D  is  a  diagonal  matrix  different  from 
the  diagonal  matrix  D.    Solving  for  C  by  matrix  methods  yields: 


C«  A"'D 


(2-10) 


For  A     to  exist,   the  determinant  of  A  must  not  be  identically  equal  to  zero. 
This  requirement  is  also  necessary  for  n  inputs  to  specify  n  outputs  independently. 
Thus,  a  condition  where  |A|=  0    implies  a  lack  of  independence  among  the 
system  equations  and  the  impossibility  of  completely  controlling  n  outputs   U'L 
by  n  inputs    01    • 

Equation  (2-10)  can  be  rewritten  as: 

Cij  =  y  &"!:*.  dfcj  (2-n) 

D  is  a  diagonal  matrix  with  elements  d..  and  therefore  the  terms  in  the  summation 
are  non-zero  only  for  k  =  j.     Thus, 

Cij  =    CUjdjj  (2-12) 

Considering  the  case  where  j  =  i,  equation  (2-11)  reduces  to  C-  -  OC^dn     or: 


da  =  -%  (2-i3) 

CCa 

Note  that  equation  (2-13)  is  the  open  loop  transfer  function  of  the  overall  system. 


Here    &  £/    is  known  and  c..  is  designed  so  that  the  closed  loop  transfer  function, 
,  provides  favorable  response  characteristics. 


After  the  c^  have  been  found  by  conventional  compensation  methods,   the 
off  diagonal  elements  of  the  controller  matrix,  c^,  are  found  by  combining 
equations  (2-12)  and  (2-13)  to  obtain:  , 

J  rrl.  i  J-  I  *»n.  J 

If  any  of  the  CCjj  =0,   then  the  choice  of  cL.  is  arbitrary.    It  follows  that 
the  off  diagonal  components,  c,   are  then  determined  from  equation  (2-12). 


i-i 


Although  d..  is  arbitrary  when  &  .;  =0  ,   it  should  be  chosen  such  that  good  closed 
loop  response  characteristics  are  obtained.    Thus,  it  should  be  chosen  fairly 
large  in  order  that  outputs  approximate  inputs  with  a  minimum  steady  state 


error*    This  may  be  seen  from  the  relation,  Lfi  =  --^-X^  •    Naturally,  the  d 


,  dii 


can  not  be  overly  large  otherwise  the  c.  must  assume  impractical  proportions 
through  equation  (2-12). 

In  summary,  the  c. .  are  determined  by  the  design  of  the  open  loop  transfer 

function,  6il  =.  ..T^     where    CC  r,     is  known.    Thus  the  c.  are  essentially  compen- 

all 
sation  so  chosen  that  is  stable  with  favorable  response  characteristics. 

The  off  diagonal  elements  of  the  controller  matrix  are  then  found  from  equation 

(2-14).     Provided  no  difficulties  are  encountered  in  compensating  the  d.,,  the 

system  performance  should  be  as  desired,  namely,  non-interaction  of  control 

variables. 
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Chapter  III 
Error  Analysis 

There  are  essentially  two  basic  types  of  system  errors.    The  first  type 
of  error  is  a  dynamic  error  resulting  from  random  disturbances  entering  the 
system.     The  second  type  are  those  errors  introduced  into  the  system  through 
the  A  and  C  matrices  due  to  inaccuracies  in  describing  the  true  physical  system 
and/or  inaccuracies  in  the  synthesis  of  the  controller  elements  or  transfer 
functions.     Errors  of  the  second  type  could  obviously  result  in  some  system 
interaction  and  conceivably  stability  problems. 
Disturbance  Errors 

In  analyzing  errors  of  the  first  type,   the  disturbances  were  considered 
as  input  signals  to  the  A  and  C  matrices  as  shown  in  Fig»  2.    It  is  supposed  that 
m  disturbance  signals,  U^ ,   are  introduced  into  the  C  matrix  and  r  disturbance 
signals,  aJa    ,  are  introduced  into  the  A  matrix.    These  input  signals  may  be 
represented  as  an  m  dimensional  vector,  La  ,  and  an  r  dimensional  vector,  V  . 
As  an  example  of  the  source  of  these  inputs,  the  U.   type  of  input  might  be 
considered  as  Tl    noise  or  drift  type  signals  introduced  at  the  inputs  of  the  YL 
c.  components.     The    V    type  of  signals  are  introduced  into  the  physical  system 
being  controlled  and  hence  could  be  torque  loads,  gust  effects,  etc.    Regardless 
of  their  source,  it  is  desirable  to  minimize  the  effect  of  the  La.  and  V  dis- 
turbances. 

It  will  be  assumed  that  LA  and  V    are  related  to  the  system  in  vector- 
matrix  form  as  follows: 
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Fig.  2    Introauction  of  Disturbances  into  the 

Control  System 
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S=  Ce  +■  c*cc 

a  -  A£  +■  A*V  (3.1} 

In  the  above  equations,   C*  is  an  n  by  m  matrix  and  A*  is  an  n  by  r  matrix. 
Noting  that  £=  X-L/ ,  the  system  error  equation  can  be  found  from  equations 
(3-1)  to  be:  * 


e=+[x+Ac]~,(x-Ac*u.-A*v) 


(3-2) 


Specifically,  the  errors  due  to  the  U.  and    V  inputs  can  be  expressed  as 

€u=   -  [r+Acl",AC*U 

,  (3-3) 

6V  =    -  [i+ac]"  A*V 


In  designing  a  control  system,  the  designer  will  have  to  investigate  the 
errors  resulting  from  the    U.    and  V   disturbances.    This  will  require  an  investi- 
gation of  the  above  vector  equations  to  determine  the  error  associated  with  the 
specific  input  disturbance.    This  may  be  illustrated  by  considering  the  effect 
of  a  disturbance  input  to  a  single  component  of  the  C  matrix  in  a  two-dimensional 
system.    Assuming  the  disturbance  as  an  input,  U,  ,  to  the  c12  component  of 
the  C  matrix,  equation  (3-3)  becomes: 


or 


v -*«<r  [t  a  ft]  »■ 
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As  developed  previously,  for  non-interaction,  AC  is  a  diagonal  matrix  which 
was  denoted  as  D.     Therefore, 


and 


e,  =  - 


i+di, 


o 

>»c,o 

1 

(a„c«) 

a. 


ti,.   = 


^  - 


i+d. 


ia 


a*,  c(J 


u., 


U.I 


14-d 


tz 


For  illustration  purposes,  c,       will  be  the  only  error  considered.    From 


equations  (2-13)  and  (2-14): 


d-ii  = 


-i 


ail 


;    Ctz  «  i£L  Cm  =  -  &*  c 


a 


-i 


a, 


2J. 


a-,  o« 


Therefore, 


e,   =  a*c 


22 


'«., 


1  +  ^/ar,', 


Assuming  the  hypothetical  system  represented  by: 

.  z 


A- 


l 
P 

.2 


pCptO 


.+  1 


then, 


P    '     fli  "     pCp*.)*       p+l 


a-',       4^     pTp+O 
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Now,  if  the  control  matrix  diagonal  elements  are  designed  as: 


C„=  i 


then, 


Cl2<  =  -.z 


r-^-i  • 

pfp+OJ     ) 


L?\  —  —>  z 


'21 


Therefore, 


t 


u'=     n 


P+IV 


«., 


Since  G-X-W,  then  the  output  y,,  with  x    assumed  as  zero,  will  be: 


«*- 


.2 


;^h 


It  is  readily  apparent  that  this  transfer  function  presents  an  undesirable 
situation,  since  at  low  frequencies  the  magnitude  of  the  output  will  be  large 
because  of  the  term  /2/p.    In  addition,   the  output  diverges  with  time  for  a 

step  input,  U,  „    Thus,  in  this  hypothetical  case  the  designer  would  be  forced 

i 

to  modify  the  controller  elements  in  order  to  provide  good  system  response 
and  at  the  same  time  minimize  the  disturbance  errors,   which  may  be  difficult, 
if  not  impossible. 

One  method  of  reducing  the  effect  of  the  disturbances  which  may  prove 
effective  in  a  practical  design  is  to  introduce  a  high  gain  prior  to  the  disturbance 
inputs,,     Fig.  3(a)  illustrates  a  basic  circuit  for  a  scalar  system  (one  of  the  non- 
interacting  subsystems  for  example)  where,  with  X-0  ,  the  transfer  function 
relating  the  output,     U     ,  to  the  disturbance  input,     LU    ,   is: 


LL 
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Fig.  3(c) 
A  Metho;i  for  Reducing  the  Effect  of 
Disturbances  Entering  the  Control  System 
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a  =  XkL-  =  Up) 

*  l+Tc(p)  (3-4) 

With  the  introduction  of  a  high  gain  element   AX     as  in  Fig,  3(b),  the  circuit 
can  be  reduced  to  that  shown  in  Figo  3(c).    From  this  latter  figure,  with  X=-0  • 

a       i+r<p)A  (3"5) 

By  proper  design,    AjC   can  introduce  high  gain  over  the  desired  frequency 
range  and  consequently  reduce  the  output  due  to  the  disturbance  inputs.  y/X. 
is,  of  course,  designed  to  maintain  good  system  characteristics  for  U   with 

X    as  an  inputo 
Component  Errors 

To  analyze  the  second  type  of  error  it  will  be  necessary  to  refer  to 
equation  (3-2)  which  is  repeated  below: 


6=    [x+ac]"1  (X-AC*a-A*V) 


(3-2) 


The  elements  of  the  A,  A*,   C  and  C*  matrices  can  be  redefined  in  terms  of  a 


lowest  common  denominator  as  follows: 


a*j  =    a' (pi 

Ctj   =     Ci(  (p) 

dc  (p)  (3-6) 

clj  =    £ii£i 
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^  ,~~.  /^t  »  O?  * 


where  O-cj  ,  Cq  ,  &u  ,  Cuj  are  polynomials  in  the  operator  p  and  da(p)  and 
dc(p)  are  the  lowest  common  denominators  of  the  a^  and  c^  components  re- 
spectively., 

From  the  above  definitions,  it  is  apparent  that  the  a*,  have  the  same  poles 
as  the  a.,  by  definition.    The  same  is  true  for  c*   and  c.o     This  generally  is  the 

case  and  if  it  is*not,  then  the  nonacceptable  poles  may  be  removed  and  used  to 

I 
define  a  new  vector,    (X     or   V    .    It  is  now  possible  to  define  new  matrices 

whose  elements  are  polynomials  in  the  operator  p.    Thus, 


ff-  (21;)      ;      ff-.W) 


(3-7) 


Note  that    A  =   4-  A  .    A*  =  -f  A  *  ,    C  *    '     c    *     C*  ^  4-  C*         Hence, 
da,        ^  da  Oc  dc 

equation  (3-2)  may  be  written  as: 

£=[cUcX  -f-AcJ'^a^X-ACV-^AN)         (3_8) 


^y       o^jt     /N'       /V» 


In  the  above  equation,  the  components  of  A  ,  A   ,  C  ,  C  ,  0^  and  Oc 

are  polynomials  in  the  operator  p„    Hence  any  poles  of  equation  (3-8)  result 

O^dc  X  "I-  A  C  „    If  determinant     dadc  I  -*-  A  C  I 

has  no  zeroes  in  the  right  half  p  plane,  or  on  the  jw  axis  (imaginary  axis),  then 

the  system  is  stable  for  disturbance  inputs  as  well  as  normal  inputs,    X^   • 

However,  if  determinant   |cLdVlI"J~AC        and  d  d    have  zeroes  which  are 

'  ■  a  c 

common,   and  if  these  zeroes  are  in  the  right  half  p  plane  then  the  system  will 
be  stable  for  X     inputs  but  could  be  unstable  for  (X  and  V  inputs.    If  dadc 
has  no  zeroes  in  the  right  half  p  plane  and  the  system  is  designed  so  that 
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.  ,  j      and    — —    are  stable,  then  the  error  vector,  c:  ,  must  be  stable  for  all 
inputs, X  i  U.  and  V  . 

When  there  are  small  changes  in  the  coefficients  of  A  and  C  due  to  inaccur- 
acies in  the  synthesis  of  C  or  inexact  knowledge  of  the  coefficients  of  the  equations 


dadcX  -hAC 


will 


of  the  physical  system,  the  zeroes  of  the  determinant  of 
change  because  t)f  the  coefficient  change.    Hence,  the  poles  of  equation  (3-8) 
will  be  altered  slightly  for  small  changes  in  the  coefficients.    Thus,  if  the 
system  is  designed  stable  and  if  dQd„  has  no  zeroes  in  the  right  half  complex 
plane  we  are  assured  that  the  system  remains  stable. 

The  small  changes  in  the  coefficients  of  A  and  C  will  result  in  interaction. 
However,  if  the  deviations  are  truly  small,  then  the  resulting  interaction  will 
be  negligible.    This  is  due  to  the  fact  that  the  d^  or  off  diagonal  elements  will 
be  small  and  hence  will  have  little  effect  on  the  overall  system.    Thus,  it  may 
be  seen  that  interaction  due  to  these  slight  deviations  can  be  minimized  only 
through  careful  attention  to  the  analytical  description  of  the  physical  system 
and  accurate  design  and  mechanization  of  the  controller. 
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Chapter  IV 
Practical  System  Problems 

In  the  previous  chapter,  the  effect  of  small  deviations,    Q    ,  in  a.. 
and  c..  were  discussed.    In  this  chapter,  the  problems  or  consequences 
attendant  with  gross  deviations  in  c^  because  of  design  advantages  will  be 
investigated.    1%  addition,  problems  arising  from  impossible  or  impractical 
compensation  requirements,   as  may  be  dictated  by  non-interaction  theory,  will 
be  discussed  and  means  of  alleviating  the  problem  pointed  out. 
Problems  arising  from  gross  deviations  in  Cj-„ 

Suppose  that  in  synthesizing  complicated  controller  transfer  functions, 
c^,  simplifications  are  made  by  neglecting  terms  in  the  numerator  and 
denominator  of  c--  that  have  "break"  frequencies  well  above  the  natural 
frequency  of  the  system.     Further  suppose  that  all  these  discarded  terms  have 
negative  roots.     The  obvious  consequence  of  such  gross  deviations  from  the 
desired  c,-  is  system  interaction  at  the  higher  frequencies.    In  addition,  one 
must  consider  the  effect  of  these  departures  on  the  stability  of  the  overall 
system. 

To  investigate  these  ideas  more  fully,  assume  a  two-dimensional  control 
system  in  which  the  exact  controller  transfer  functions,  c12  an<^  c2l  are 
altered  by  a  factor  q  to  produce  new  relations  c'     and  c'    .    Analytically, 

1   —  Z  1 


C?  =  4,  c 


f,    C'^  (4-1) 


C'2I    .  fC 


(4-2) 
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For  the  theoretical  non- interacting  system,  the  following  equations  have  been 
shown  to  exist: 

du  =   a„  C,   -h  Q,z  Czi    a  Yot  -     C"/d-\  (4-3) 

du  =  a2l  Cz  -hdvCv  a_Kx  -  C2^-a'z  (4-4) 

«/i  -    (2//   Cz    f  fl/i  C22.     =  0  (4-5) 

^2/    -  dZ\  C„   -h  Ciz2.  Cii    =  0  (4-6) 

Substituting  equations  (4-1)  and  (4-2)  into  equations  (4-3)  through  (4-6)  where 
applicable  and  manipulating  and  rearranging  where  desired,  the  results  are: 

du  =  d„  c,  +  aiz  fc  Cz,  =  />,  +  f*a*  a%-<]c<  <4-7> 

d'x  ^duCzz-h  ty  CUt  Cz  ^[a^  +  fr  a*>  a"'ytii]czi.  (4-8) 

d,'z  -  dn  fy  Ol  *hd*  Czz  =  [CL,Z  (l-p  )  CziJ  j=  0      (4_9) 

d'zi  =  An  C,  +  0%  Czi  dzz  =  \_Ciu  (i-p)  C,J  ^  0      (4-i0) 

where  the  d'    are  elements  of  a  new  and  non-diagonal  matrix,  D' . 
ij 

The  interaction  is  evident  from  equations  (4-9)  and  (4-10).    The  changes 
in  the  transfer  functions  d-  due  to  the  factors  q^  may  be  determined  from 
equations  (4-7)  and  (4-8). 

In  investigating  the  effects  of  gross  deviations  on  system  stability,  it  is 
best  to  look  at  the  system  equation  in  vector-matrix  notation.    It  is: 
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1 


r    [ifD'fD'x 


where  D'  is  a  non-diagonal  matrix. 

The  stability  of  the  closed  loop  system  is  determined  from  the  determinant 
of  |I+D  I  which  is  the  characteristic  equation  of  the  overall  system.    It  can  be 
shown  that: 

|x-KD'|  -  O+^Xi+JttJ-d/idi,  (4-ii) 

It  follows  that  if  this  resulting  polynomial  has  zeroes  in  the  right  half  complex 
plane  the  overall  system  is  closed  loop  unstable. 

One  means  of  determining  the  stability  of  the  closed  loop  system  without 
solving  for  the  roots  of  the  polynomial  is  as  follows: 


(/  +  d,'i)0+d«)-  dlz6'u    =  o 


or  when  rewritten: 

i-   EtaiL 


0 +<*;,)(' +<*») 


1L 


O+diiKi+djJ 


=  o 


(4-12) 


Stability'  is  assured  if  none  of  the  above  factors  have  zeroes  in  the  right 
half  p  plane. 

The  open  loop  transfer  functions,  d'     and  d'    ,  which  are  given  in 
equations  (4-7)  and  (4-8),   should  be  plotted  in  Nyquist  plots  and  Nyquist  stability 
criteria,  N  =  P  -  Z  applied  to  determine  if  there  are  zeroes  of  the  1  -f  d^   in 
the  right  half  p  plane. 
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Assuming  that  there  are  no  zeroes  of  the  1  +  d! .  in  the  right  half  p 

plane  and  that  d'     and  dk,  have  no  poles  in  the  right  half  p  plane,  then  for 

du  0/21 


system  stability,   the  absolute  value  of 


must  be  less  than 


unity  for  all  frequencies,  CO    .     However,   reference  to  the  Nyquist  plot  of 
a  typical  d*--  shown  in  Fig.  4  provides  the  following  important  relationship: 


i+da  £  a 


where 


a= 


ML  + 


(4-13) 


It  follows  from  equation  (4-13)  that: 

Thus,  if  (1  +  M1)  (1  +  M  )  d^2  d^  <    1  for  all  frequencies,  UJ  ,  and  (H-  d^) 
as  well  as  (1  +■  d'    )  have  no  zeroes  in  the  right  half  p  plane,  then  the  overall 
control  system  is  closed  loop  stable. 
Compensation  Problems 

It  is  conceivable  that  the  co- factors  of  the  a.,  of  the  physical  system 


in 


may  have  zeroes  in  the  right  half  p  plane  so  that  the         :      will  have  poles  i 
the  right  half  p  plane.    From  the  basic  open  loop  non-interaction  relation, 
Oft  —  — —  ,  it  follows  that  it  may  be  impossible  or  impractical  to  compensate 

da 

with  c.,.  so  that    is  stable.    In  this  event,   it  would  be  impossible 


CC- 


ii 


a 


H-dtf 


to  obtain  a  practical  non- interacting  control  system  using  the  design  techniques 
discussed  in  the  theory  section.  As  an  example,  consider  the  following  simple 
system: 


A  = 


p 


■H 


.& 

p+l 

-.OIp*l 

p-M 

•.  i'(i) 
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radium 


M  circle  _/ 
for  M  =  A-//  /^j 


dji  p/c 


one 


a  =  //?e   c//s  fence    -from   the 
critical  pomt  to  the  in- 
tersection,  of  the  M  circle 
and  the  real  axis. 

^~ 

/?«?/     <7X/S 


/ocus  ot  d-i 


Fig.  4    Nyquist  Plot  of  a  Typical  Open  Loop  Transfer 

Function,  cL  • 
11 
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Then: 


oi 


-I       __     (-.0/  p-hl)( p-hl) 


u 


-.Olp  -h.  9G> 


and: 


a 


-•  Olp  -h.  9U 


"n         i-oip-hi)(p+i) 


The  locus  of  — —  for  p  =  jw  is: 


VO-0 


But  for  stability,  N  =  P  -  Z  and  in  this  case  P  =  1»  Therefore,  N  must  be  f  1 
if  Z  is  to  be  zero.  However,  from  the  general  non-interaction  theory,  d,,=  •=^- 
It  is  quite  obvious  that  the  compensation  of    —7    such  that  N  =  +•  1  and  Z  =  0  is 


cc 


II 


a  difficult  problem.    One  means  of  compensation  is  to  use  the  prediction  oper- 

p 

ator,  £     .    This  will  provide  the  necessary  encirclements  but  is  entirely 

impractical  as  a  means  of  compensation  in  a  control  problem.    Other  than  this, 
no  compensation  could  be  found.    One  might  consider  mechanizing  the  compensa- 
tion so  that  it  cancelled  the  objectionable  poles  and  zeroes,    However,  this  wiU 
not  work  because  the  poles  and  zeroes  of  the  physical  system  are  not  known 
exactly  and  unless  we  have  exact  cancellation  the  poles  and  zeroes  are  not 
really  cancelled.    Thus,  it  appears  that  the  straight  forward  design  techniques 
developed  in  the  theory  can  not  be  used  when  confronted  with  such  a  compen- 
sation problem. 
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In  some  cases  the  foregoing  design  problem  can  be  avoided  by  merely 
interchanging  the  rows  of  the   o    vector.    This  implies  an  interchange  of 
controls o    If  the  design  problem  is  removed  by  such  a  change,  one  may  continue 
to  design  the  control  system  in  accordance  with  the  general  non-interaction 
theory. 

As  an  exarfiple  of  the  foregoing,  consider  the  same  system  of  the  previous 
example  but  with  the  o  vector  rows  interchanged  to  form  a  vector  o  and  a 
corresponding  change  of  the  columns  of  A  to  form  a  new  matrix  A'.    Thus, 


A  = 


.  z 

i 

p+l 

p+l 

■Olp-H 

.  2 

p-H       p+l_ 


;  '-( 


The  new  d. .  are  now: 


d„  =  ro.  -  cdAi  . 


a 


'-i 


p-M 


<U=X  -  £dfl  = 


a 


-  Czz.  (.2)(-  0<p-h.94>) 
p-hl 


There  are  no  poles  in  the  right  half  p  plane  of  the  above  open  loop  transfer 
functions.    Hence,   P  =  0„    In  addition,  the  design  of  cl;,and  c22  so  as  to  have 
N  =  Z  =    0  in  the  closed  loop  system  is  not  anticipated  to  present  any  real 
compensation  difficulties.    Thus,  the  impossible  or  impractical  compensation 
problem  has  been  avoided.    It  should  be  emphasized  that  interchanging  the  rows 
of  the    $     vector  will  not  alleviate  this  particular  compensation  problem  in  all 
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instances.    This  is  obvious  from  the  foregoing  example  if  either   ~TJ     or 

6CI2 

I 
"TT    have  one  or  more  poles  in  the  right  half  p  plane.    Naturally,   the  presence 
"■zi 

of  such  poles  does  not  imply  that  the  can  not  be  practically  compensated 

so  that  Z  =  0  but  rather  that  difficulties  in  compensation  are  more  probable 
than  the  instances  when  P  ■  0. 

In  a  more  Complicated  system  where  there  are  n  control  inputs  and  out- 
puts and  n  is  greater  than  two,  the  poles  of result  from  the  co-factors  of 

the  A  matrix.    These  complicated  co-factors  tend  to  increase  the  probability 
that  — r*   will  have  one  or  more  poles  in  the  right  half  p  plane.     Thus,  one 
might  expect  some  difficult  and  impossible  compensation  problems  when  applying 
the  general  non-interaction  theory  to  complicated  physical  systems „    If  an  inter- 
change of  the  controls,  (interchange  d'L )  does  not  avoid  these  problems  then  a 
different  design  approach  must  be  taken,, 

The  theory  behind  one  approach  and  the  effect  on  system  performance  is 
most  easily  described  for  a  simple  two-dimensional  system.    For  instance, 
suppose  that  &,.,  can  not  be  properly  compensated  by  c,,.    The  result  is  that 
the  system  can  not  be  designed  entirely  within  the  non-interaction  theory.    How- 
ever, a  practical  system  can  be  found  if  the  designer  will  accept  some  inter- 
action.    Furthermore,   some  of  the  non-interaction  techniques  can  be  used  once 
the  difficulty  has  been  removed. 

In  finding  a  practical  system  it  is  probably  best  to  start  with  the  basic 
expression  for  d,,«    This  is: 

<f«  =  0,.  0,  +  6UG,  <4"14> 
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In  as  much  as  it  is  impossible  or  impractical  to  compensate  cL,  with  c, , , 
then  c.,  may  in  general  be  chosen  arbitrarily.    However,  in  some  instances  c21 
may  be  chosen  completely  arbitrary  and  hence  the  arbitrariness  of  c, ,  will  be 
restrained  to  a  certain  type  of  transfer  function  which  results  in  1  +  d'     having 

no  zeroes  in  the  right  half  p  plane    In  addition,  certain  simplifying  assumptions 

I 
in   ~tt   may  result  in  restrictions  one...    In  any  event,  c..  assumes  a  new 

description  asc'    .     Since  it  may  be  necessary  to  change  c21  from  the  transfer 

function  dictated  by  non-interaction  theory,  it  now  becomes  cl,.    Therefore, 

I  f  d,',    =     I  +  dlt  Cm   4-  d,z  Cu  (4_15) 

The  procedure  is  to  make  lfd'     so  that  there  are  no  zeroes  in  the  right 
half  p  plane.    Note  that  1  +  d^  will  have  zeroes  in  the  right  half  p  plane  because 
^11  *       /OC^         can  not  ke  properly  compensated  so  that  Z  =  0. 

In  as  much  as  the  a.,  are  fixed  by  the  physical  system  it  is  much  easier 
to  design  c'     and/or  c*     so  that  there  are  no  right  half  p  plane  zeroes  in  the 
resulting  1  +■  d'   .    As  a  first  approximation,  write  c'     in  the  form  dictated  by 

XX  Z  X 

non-interaction  theory,  which  is: 

C*    -  —^  U  (4-16, 

In  as  much  as  c'     is  arbitrary  in  some  sense,  make  it  as  simple  as  desired. 
Then  substitute  into  equation  (4-15)  and  apply  Routh's  criteria  to  the  resulting 
polynomial  to  determine  whether  or  not  these  are  zeroes  in  the  right  half  p  plane 
If  this  test  fails,   alter  cl.  by  a  factor,  a,  where  a  is  a  polynomial  in  the  operator 
p  or  a  ratio  of  polynomials,,     Again,  test  the  resulting  polynomial  for  positive 
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roots.    Eventually,  a  q  will  be  chosen  such  that  1  -+■  d'     has  no  zeroes  in  the 
right  half  p  plane.    In  addition,  if  q  is  chosen  so  that  the  "break"  frequencies  are 
above  the  system  natural  frequencies  then  the  interaction  will  be  restricted  to 
these  higher  frequencies. 

Having  removed  the  difficulty  in  d     ,   attention  can  be  focused  on  d2„.    An 
attempt  should  be*  made  to  design  it  in  accordance  with  the  general  theory  of  non- 
interacting  systemSo    If  no  compensation  difficulties  are  experienced  then  c 
and  c12  will  be  dictated  by  the  general  theory.    In  addition,  1  +  d      r  1  +•  dl? 
and  there  will  be  no  right  half  p  plane  zeroes.    If  compensation  difficulties  are 
encountered  then  follow  the  procedure  outlined  for  d'  J0 

Naturally,  the  system  has  some  interaction  which  appears  through  the 
d'     and  d'     transfer  functions.    In  fact,  for  the  case  presented  here 

d>,'2  =  a„  cU  +  a,*  cL  =  o   [^c«  ^*  J*j        (4-i7) 

di  =  a2lcl,  +  &x  c'zl  ?o  (4-i8) 

Due  to  interaction,  the  overall  system  closed  loop  stability  is  dictated 
by  the  zeroes  of: 

(  I  i-d'n  )  (  I  +  d£i  )  -    d,2  d*i      =   O  (4-19) 

Determining  the  stability  through  this  equation  has  already  been  discussed  at 
length  in  the  previous  section.    In  the  event  that  there  were  zeroes  of  equation 
(4-19)  in  the  right  half  p  plane  they  would  stem  from  the  interaction  terms, 
d'12  and  d^..    Thus,  it  would  be  necessary  to  choose  another  q  factor  so  that 
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both  the  1  +•  d'. .  and  I-0'2-""      contributed  no  right  half  p  plane  zeroes  to 

11  (/+d,',)(i+d»'*) 

equation  (4-19). 

In  the  case  considered  here  di2.— O  .    Thus,  system  stability  is  assured 
if    (  l  4-  On  )  \  I  +"  diz.  )  have  no  zeroes  in  the  right  half  p  plane.    Actually, 

Qiz  =  0       due  to  the  slight  deviations  in  a.,  and  c.  of  equation  (4-17).    How- 
ever,  it  is  not  expected  that  such  slight  interaction  through  this  term  would  cause 
instability  when  coupled  with  d' 

In  summary,  when  compensation  difficulties  arise  and  interchanging  the 
rows  of  the    o     vector  does  not  alleviate  the  problem,   then  it  is  necessary  to 
deviate  from  the  non-interaction  theory  to  the  extent  that  there  are  no  zeroes 
in  the  right  half  p  plane  of  the  resulting  characteristic  equation  of  the  system, 
namely: 


(\+d'u)0+dU)  -d!idi    ^o 


Naturally,  interaction  will  be  introduced  through  the  d'..  transfer  functions.    In 
the  case  considered  here,  the  amount  of  interaction  as  a  function  of  frequency 
can  be  found  from  the  following  transfer  functions: 

</»/*  =   ik  (4"21) 

1  /X'  |D'I 


where 


l^'l  »  6 +<&)('+*&)-<&<& 
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Chapter  V 
The  Synthesis  of  an  Aircraft  Control  System 

In  present  day  aircraft,  the  auto  pilot  performs  many  functions  in  the 
integrated  control  system.    Although  the  system  inputs  are  generally  restricted 
to  combinations  of  elevator,  rudder,  aileron  and  throttle,  the  system  controlled 
outputs  are  many?    For  instance,  an  auto  pilot  system  could  control  the  aircraft 
velocity  vector  through  the  control  of  the  pitch,  yaw  and  roll  rates.    Other 
systems  could  control  either  altitude,  angle  of  attack,  airspeed,  glide  path, 
Mach  number,  lateral  acceleration,  turning  rate,  etc. ,  or  combinations  of  the 
foregoing.    It  is  the  intent  here  only  to  point  out  that  the  auto  pilot  can  be  used 
to  control  many  variables  and  each  subsystem  of  the  auto  pilot  can  control  as 
many  outputs  as  it  has  inputs. 

In  order  to  point  up  the  theoretical  analysis  and  the  practical  problems 
associated  with  a  design  based  on  the  general  theory  of  non- inter  acting  controls, 
an  auto  pilot   subsystem  shall  be  synthesized  in  this  section.,    Naturally,  the 
goal  of  such  a  study  is  to  arrive  at  a  controller  design  that  permits  good  system 
performance  as  well  as  the  desired  non-interaction.    However,  as  was  pointed 
out  in  the  preceding  sections,  a  perfect  non-interaction  over  all  frequencies  is 
not  obtainable  nor  is  it  necessarily  requiredo    In  addition,  complex  physical 
systems  such  as  an  aircraft  could  easily  require  impractical  or  impossible 
controller  compensation  when  applying  the  general  non-interaction  theory. 
Thus,  the  designer  may  have  to  accept  some  interaction  in  order  to  preserve 
system  stability  o     Furthermore,  he  may  be  forced,   at  least  in  part,  to  deviate 
from  the  general  theory  in  order  to  design  a  practically  realizable  system  with 


31 


desirable  characteristics. 

The  subsystem  to  be  synthesized  is  one  dealing  with  longitudinal  control 
in  which  the  inputs  are  elevator  and  throttle  deflections.     The  controlled  out- 
puts are  glide  path  angle,   Qw   ,  and  indicated  airspeed,  /"U7>    ,     It  should  be 
emphasized  that  this  study  was  chosen  only  as  a  means  of  lending  practical 
significance  to  the  theory  and  is  not  considered  as  a  demonstration  of  a  specific 
auto  pilot  design. 

The  assumed  aircraft  and  flight  condition  is  an  F-86   'Saberjet'  making 
an  approach  down  a  3°  glide  slope.    Small  perturbations  from  symmetric 
equilibrium  flight  will  be  assumed  thus  permitting  the  use  of  the  linearized 
longitudinal  equations  in  the  ensuing  analysis.    Upon  determination  of  the 
controller  design  for  the  best  system  performance,  the  synthesized  system 
will  be  studied  on  the  differential  analyzer.    Such  a  study  should  clearly  demon- 
strates the  applicability  of  the  theory,  error  analysis  and  the  practical  system 
problemSo 

The  large  motion  longitudinal  equations  of  motion  for  an  aircraft  using 
a  combined  body- wind  axis  system  are  shown  in  Ref.  5  to  be: 


•  a  # 

f^x=0  =  ->nVp  -»-Fwx  =  -wVp-t-hxCoso'  -*-rs-3irW  (5-1) 

Fwa  =  o  «  ~™VpQ»i-Fv<*  =  yY)VpQ»-Pxsir)o(  +  p2CQ$ot  (5-2) 

4-  2.z  +-  ma  cos  0w 


(5-3) 
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The  linearization  of  these  equations  is  carried  out  in  Appendix  I.    In  addition, 
a  fourth  linearized  equation  is  obtained  so  that  there  are  four  simultaneous 
equations  involving  the  four  output  variables,  AJp  ,<=>/,   #•  ,  6^  and  the  inputs, 
dj  and  de  . 

It  is  desired  to  determine  the  transfer  functions  that  relate  /vJ^  and  d^to 
Sf  and  £g  „    It  stlould  be  clear  from  the  theory  that  these  transfer  functions 
are  the  a.,  elements  of  the  A  (aircraft)  matrix.     The  mathematical  steps  taken 
to  determine  these  elements  from  the  four  simultaneous  equations  are  completely 
outlined  in  Appendix  I. 

The  elements,  a^,   are  merely  a  ratio  of  polynomials  in  the  operator  p 
or  Do    However,   they  derive  their  existence  from  the  coefficients  of  the  line- 
arized longitudinal  equations „     These  coefficients  are  composed  of  thrust  terms, 
gravity  terms  and  aerodynamic  force  and  moment  terms  expressed  through  the 
aircraft  stability  derivatives  and  trim  lift  and  drag  coefficients. 

The  stability  derivatives  for  the  F-86  were  obtained  from  Ref.  6  for 
flight  at  sea  level  at  163  knots.     They  are  listed  in  Appendix  II  along  with 
other  physical  data  pertaining  to  the  F-86.     Calculation  of  the  trim  lift  and 
drag  coefficients  as  well  as  the  b.    and  f     coefficients  of  equations  (1-9) 
through  (1-12)  are  also  found  in  Appendix  II.    Knowing  the  b-  and  f«  co- 
efficients permitted  the  calculation  of  the  a-  elements  from  equations  (1-18) 
through  (1-21)  of  Appendix  I.     The  results  are: 

7l>-.Ol£)(tf>l.4a£P  +3.4771) 


a..,  =  \o.zz 


(D*  4- 1. 41 443)  4-3. 4S,34)(DN-0lSkD  4-. 0245) 
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CL.2.  =  - .  ass- 


(JP-iScQCd-H.i) 

(pz  +  l.4l44J>  f  3.4534  )(DV.0l5fcD  40*45') 


61  9,  =    .  QOfoZS 


21 


('P+-  1.409)  CdV.701:D  -t- 3.087) 


(Vf  I.4I44P  4- 3.4534) (dV.O/S&D  +.0245) 


CUt  «  .07 


(;D-.Q07J7)(l>+7.fcSXP-8'S) 

(D2+l,4i44l>-l-S.4534)(DV.OI5-bl)  K0245) 


The  two  quadratic  factors  in  the  denominator  of  the  "a". .  terms  provide 
information  as  to  the  period  and  damping  of  the  short  period  and  phugoid  motions. 
The  period  and  damping  factors  obtained  are: 

Short  period  motion  Phugoid  motion 


11  =    —•  =     3.38  sec. 
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T  =  2*    - 


LL>fM 


40  sec. 


rP=  • 


0842 


It  may  be  seen  from  Figs„  1  and  5  that  in  vector- matrix  notations: 
-  TASS 


*J 


where  —     l^Vi) 

5  =  (sa) 

A    ^    {OLLj) 
For  the  specific  problem  under  study: 

y,  s  /Op    ;  y,  =  ^   ;    5  =  ST     ;    S,  =  £T 

Therefore: 

/Op  =  l-tn  2,1  3,1 )  ST  +(+h  S,2  522)  5e  *  a„  <5r  f  a,2  &         (5-4) 
Gw  -  I  tz  #2i  S,i)  dj  4-  (  fu  #22  522 )  oe  =  ^21  or  +~  Uzz.de         (5-5) 


To  determine  the  a.-  terms  of  the  above  equations,  it  is  necessary  to  specify 
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the  elements  of  the  diagonal  servo  and  transducer  matrices,   S  and  T  respectively. 
This  has  been  done  in  Appendix  HI.    In  addition,  reasons  governing  the  specifi- 
cation of  these  elements  or  transfer  functions  are  included  in  the  appendix.     With 
this  knowledge,   the  a.,  terms  of  equations  (5-4)  and  (5-5)  are  found  to  be: 


&„   =   10.22 


a*  -  -.  8&5 


(D-.0i2.Ud2'4J.42Z:D  4-3.4771) 


#21  =..OObZS 


Q-ll-  .01 


£p2  +- 1. 4I44P  4-3.4534 )(DN-.0l5ifeD  h.0Z46~)  (2D4-I)  (•  2D-HJ 

~ (S-lSb)  (D4-(,|) 

toa  4-1. 4I44D  4-3.4554)  (DVoiSbJ)  -h,0245)C.04]>+-j)(-  2-D4  l)_ 

(Dfl.4C^)(DZ4>,7QIJ)-f-3.oe>7) 


(02- 4-|. 4l44D4-3.4S3^)(t)V.0i5bDf. 0245)  C^D4-|) 


(D--.C07i7)(D-t-7>^5)^-e^5) 

(DSl.4l441>  4  3. 4534)  (  dV  QI5feD  4r.QZ4sU-  Q^-^-H  ) 
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It  has  been  pointed  out  that  d. .  =  c    /a7.    is  an  open  loop  transfer  function 

n       ii     11 

for  one  "channel"  of  the  auto  pilot  system.     From  Fig.  5,  it  may  be  seen  that 
the  open  loop  vector-matrix  equation  is    U  =  ACG  .    It  follows  that; 


*,-YH.^-   Cv/ff,t>   C^al,|» 


(5-6) 
(5-7) 


9<A> 


It  is  evident  that  equations  (5-6)  and  (5-7)  are  the  open  loop  transfer  functions 

for  the  velocity  and  pitch  angle  "channels".    It  is  desired  that  these  channels 

be  closed  loop  stable  as  well  as  having  desirable  response  characteristics.    An 

attempt  is  made  to  design  the  controller  transfer  functions  c      and  c      so  that 

&  11  22 

these  requirements  are  met.    If  they  are  satisfied,  then  the  other  controller 
transfer  functions,  c]2  and  c„    are  determined  from  equation  (2-14). 
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Fig.  5  The  Open  Loop  Control  System  Configuration 
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The  determinant  of  the  A  matrix  was  found  to  be: 


A 
A 


alla22      a12a21 


=  .71 


(D^7.74XD-8.S87XD^-t-l.4l5Df3.4Sfe)(D%.Ql87])-»->Q£l^>) 


Substituting  into  equations  (5-6)  and  (5-7)  and  cancelling  equal  or 
essentially  equal  terms  in  numerator  and  denominator  (excepting  terms  in 
right  half   plane  )  provided: 


Y0,  =  C„  Oo.z) 


(p-6.S£7) 


(p-.oo7»7)(p-8,s)fepn)C2p+-0 


1o4=Ctt  (1.333) 


(.ia<?p+l)  (.nfcSp-Q 

(p-.0«2)(.  298  p*>.  4<X?p  *-|  ) (.04pf  I  ) 


(5-8) 


(5-9) 


In  deciding  what  c,  ,  should  be  in  order  that   ^ —     be  closed  loop  stable 

with  good  response  characteristics,  Nyquist  plots  had  to  be  used.    A  sketch 

V 
of  the  Nyquist  plot  of      °i      is  shown  in  Fig„  6.    There  are  two  poles  in  the  right 

Cnk 


half  p  plane  of 


Nyquist  stability  criteria  indicates  that  N  =  +■  2  for 


Z  =  Oo    This  implies  that  c...  must  be  chosen  such  that  the  Nyquist  plot  of  the 
resulting  Y      or  d^  must  encircle  the  critical  point,  -1,  0,  twice  in  a  counter 
clockwise  direction.     Without  doubt,  the  required  compensation  presents  a 
rather  formidable  problem.    In  fact  no  practical  compensation  could  be  found. 
Thus,  it  is  clear  that  the  system  will  be  closed  loop  unstable  in  that  it  can  not 
be  properly  compensated  by  c^. 

In  an  attempt  to  alleviate  the  above  design  problem,  the  rows  of  the  o 
vector  were  interchanged.    As  a  result,  it  was  possible  to  compensate  d      in 


accordance  with  the  general  non-interaction  theory.    However,  compensation 
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Fig.  b    A  Sketch  of  the  Nyquist  Plot  for  the  Transfer 

Function,  X/C„  K 


38 


problems  arose  in  Y      or  cL2  and  no  practical  means  of  compensating  were 
found.     Thus,  this  method  failed. 

The  next  approach  was  to  alter  d..  with  a  factor  q  such  that  there  were  no 
zeroes  of  1  +  d'  .  in  the  right  half  p  plane.    This  approach  is  discussed  at  length 
in  the  section  on  system  problems.    From  that  section,  it  is  known  that  a  more 
general  form  for  1  +  Y      or  1  +  d'     is  given  by  equation  (4-15)  which  is  repeated 
here: 

l-l-d',   -  Aid   +  fliiCsi    -hi  (4-15) 

It  was  pointed  out  in  Chapter  IV  that  c\,  or  ck,  could  be  somewhat 
arbitrary  when  d      could  not  be  properly  compensated  by  the  non-interaction 
design  methods.    In  the  case  at  hand,  instead  of  making  c*     completely  arbi- 
trary, let  us  assume  that  the  p-8.  5  pole  is  cancelled  by  the  p-8.  587  zero  in 
equation  (5-8).    Of  course,  it  is  understood  that  this  is  incorrect  in  that  one 
of  the  right  half  p  plane  poles  are  being  neglected.    Nevertheless,  compensating 

the  remaining  Y^    with  P  ■  1  results  in  c'     being: 
°1  11 

Naturally,  this  is  not  the  necessary  compensation  for  closed  loop  stability  of 

—    as  dictated  by  non-interaction  theory. 

As  a  first  approximation  to  ci  ,  let  it  be  dictated  by  the  equation: 


c£,  =  -  p.  cl, 

Substituting  for  a2p  a      and  c'      results  in: 
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0  ip-hi)  Cp'-ooin)  (.i3oqp+i)Cni&p-i)        (5_11) 

Substituting  equations  (5-10)  and  (5-11)  as  well  as  a^  and  a^  into  equation 
(4-15)  and  applying  Routh's  stability  criteria  to  the  resulting  ninth  order  poly- 
nomial  indicates  that  there  are  zeroes  in  the  right  half  p  planeo     However,  by 
altering  cU,  above  by  a  factor  q2,  where  q2  =  -(.  1178p  -  ly(.  04p  4-  1),  the 
new  c21  becomes: 

s>t    _        > 002.92.  (,li(ptl)CZ24p7'+-Z2-'7p-hi) (5_12) 

ai  (p~.oonn)(.i p+i)  (,i3QCip4-i) 

Going  through  the  same  substitution  procedure  and  applying  Roth's  criteria 
indicates  that  there  are  no  zeroes  of  1  -f-d^  in  the  right  half  p  plane.    It  follows 
that  the  velocity  "channel,"  Y,,  is  now  closed  loop  stable.    In  addition,  inter- 
action is  present  in  the  system  because  the  c'     and  ci.  of  equations  (5-10)  and 
(5-12)  are  different  from  the  c, ,  and  c2,  that  would  have  been  dictated  by 
following  the  non-interaction  theory. 

The  d22  transfer  function  was  determined  from  the  non-interaction  theory. 
It  is  shown  in  equation  (5-9),  A  Nyquist  plot  of  -Y     /C22K   is  shown  in  Fig.  7. 

cL2(or  YQ  )  was  properly  compensated  so  that 


Lg£jH 


By  making  C      -T5 

22  |1c^.p+ 

there  were  no  zeroes  in  the  right  half  p  plane.    Furthermore,  the  compensation 
was  chosen  so  as  to  give  good  transient  response,, 

Having  determined  c      in  accordance  with  the  non-interaction  theory,  then 
it  was  permissible  to  find  c,2  from  equation  (2-14).    Therefore: 


2  10° 
150° 
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CIZ  =  ~4ukCu  =    -2..o+s(i-zp+i)(.<iip+i)(zpH)(.oot>7p-,)         (5_13) 

O.,,  (p-.OI2)(,Q4pt-l)C0Zpi-l)(.2tfp1--t-.40<ip  +  i) 

In  as  much  as  terms  such  as  („  00667p  -  1),   (.  02p  ■+- 1)  and  (.  04p  +  1)  are 
only  effective  at  the  higher  frequencies,  c      can  be  simplified  by  neglecting 
them.    Naturally  this  will  give  rise  to  interaction  at  higher  frequencies.    But 
this  is  of  no  real  consequence  provided  that  these  gross  deviations  do  not  im- 
pair the  system  stability  through  the  d'    d'       term  of  the  system  characteristic 

12     21 

equation,,    In  the  problem  considered  here  the  stability  is  not  impaired  by  the 
interaction  terms. 

In  light  of  the  foregoing: 

C'   «  $ /),  s    *o*s (tep+i)C*ip+i)(?.pi-i) (5_14) 

^       V  (p-.o<a)C.2Z?p*-+.*>9p+j) 

where  .  v  ,  s 

q  _     (.Q4pi-/)(.0Zp+i) 

V  -  (.oobip-  i) 

Attenuation  plots  of  c'  , ,  cnn  and  the  c     of  equations  (5-10)  through  (5-14) 

1 1       2i£  ij 

are  shown  in  Fig.  8„    One  may  readily  see  that  the  effect  of  neglecting  the  high 
frequency  terms  is  to  level  off  the  gain  at  high  frequencies,  such  frequencies 
being  above  the  natural  frequency  of  the  overall  system. 

The  amount  of  interaction  contributed  through  d'     and  d'     can  be  found 
from  equations  (4-20)  and  (4-21),    In  the  case  considered  here,  d"     and  d' 


are  as  follows: 


dt'z  =     &t  p  G*    +■   GLi  Czz.  (5-15) 


qp  ~  ja/?j/A/5  vw 
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d'  =  cui  c«  -h  cL  a 


'<?/     ~     U*l     II  z/   LX2i  (5-16) 

where  c'     is  given  by  equation  (5-12)  and  c     q    is  given  by  equation  (5-14). 
Computer  Simulation 

With  the  c..  designed,  the  auto  pilot  system  was  mechanized  for  study  on 
the  Reeves  Electronic  Analog  Computer  (REAC).     The  mechanization  used  is 
shown  in  Fig.  9=, 

System  response  to  step  inputs  of  command  velocity,  — Bu.  ,  and  command 
glide  path  angle,   #w.   ,  are  shown  in  Figs.   10  and  11.    Note  that  the  interaction 
is  essentially  zero  except  at  high  frequency.    However,  this  was  expected  in 
that  high  frequency  terms  in  c'     and  c'     were  neglected.    There  is,  although 
not  evident  in  Figs.   10  and  11,   a  low  amplitude  damped  oscillation  that  affects 
the  steady  state  character  of  — £-  and  (L. .     These  oscillations  are  the  direct 
result  of  disturbance  inputs,  (i^  ,  that  enter  the  system  through  the  c...    In 
fact,   Figs.  12(a)  and  12(b)  show  these  oscillations  in  ixEand  9^  when  there 
are  no  inputs  other  than  disturbances  in  the  system. 

As  a  means  of  verifying  the  magnitude,  frequency  and  damping  of  the  low 
amplitude  oscillation,   a  step  input,  x2,  was  introduced  into  c21  through  the  in- 
put  to  amplifier  R-8  of  Fig.  9-     The  effect  of  the  disturbance  input  on  — -*-—  and 
0^  is  shown  in  Figs.   13(a)  and  13(b)o 

An  analytical  study  of  the  system  was  made  in  order  to  determine  the 
relation  between  -£-  and  x  and  also  to  verify  the  traces  seen  in  Fig.  13  as 
those  that  stem  from  disturbance  inputs.    In  the  case  where  there  are  no  desired 
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perturbations  from  steady  state  condition,  then  s^p^o.    Therefore,  it  may 
be  shown  that: 

Substituting  into  equation  (5-17)  and  simplifying  results  in: 


[p%4?.Jp8  +<gOOp7  +  S(l70p<'  f  Z70SDp5->-c?0|4t>Op4tl3lfc)00pi  +  |7/,0UOp; 

+■  5  /  50p+-  4 1 OOJ 


where  the  static  gain  is  unity- 

Application  of  Routh's  criteria  to  the  ninth  order  polynomial  in  the  denomi- 
nator of  equation  (5-18)  indicates  that  there  are  no  positive  roots.    Therefore, 

the  system  is  stable  for  the  disturbance  input,  x  .    In  addition,  it  may  be  shown 

2 
that  (p  +•   „  0194p  +•  o  0210)  is  a  factor  of  the  ninth  order  polynomial.    Further- 
more, it  is  the  only  factor  which  has  a  low  characteristic  frequency.    In  fact 
this  frequency  is  .  0231  cycles  per  second  and  the  damping  factor  is  f  =  .O(o7  . 
It  is  quite  apparent  that  this  quadratic  factor  alone  accounts  for  the  lightly 
damped  high  period  oscillation  seen  in  Fig.   13.    Naturally,  similar  analysis 
can  be  made  for  other  disturbance  inputs  into  the  system. 

The  low  amplitude  oscillation  discussed  above  as  well  as  other  responses 
to  disturbance  inputs  can  have  a  detrimental  effect  on  the  system  due  primarily 
to  their  magnitude  and  damping.    The  addition  of  a  high  gain  loop  would  introduce 
gain  prior  to  the  introduction  of  noise  and  hence  make  the  system  less  sensitive 
to  the  disturbance  inputs. 
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Chapter  VI 
Conclusions 
lo        In  theory,  control  systems  based  on  non-interaction  of  control  variables 
are  feasible  but  somewhat  complex.    For  complete  non-interaction,  the  theory 
requires  precise  mathematical  knowledge  of  the  behavior  of  the  system  to  be 
controlled  and  exact  mechanization  of  the  controller  transfer  functions.    In 
reality,  these  requirements  generally  can  not  be  met  but  can  be  approached 
to  the  point  that  the  resulting  interaction  is  negligible  for  all  frequencies. 

2.  Small  deviations,    €'L   ,  ina..  due  to  insufficient  knowledge  of  the  physical 
system  and/or  in  c.  due  to  the  inability  to  mechanize  transfer  functions  exactly 
should  not  impair  the  stability  of  the  system  if  the  d..,  or  open  loop  transfer 
functions,   are  designed  with  a  reasonable  amount  of  stability, 

3.  Gross  deviations  in  c    ,  so  as  to  obtain  a  simpler  mechanization  can 

ij 

result  in  considerable  interaction  and  possible  instability  in  the  closed  loop 
system.    If  deviations  are  restricted  to  terms  effective  only  at  frequencies 
well  above  the  natural  frequency  of  the  system,  then  the  interaction  will  be 
restricted  to  these  frequencies.    However,   system  stability  can  be  impaired 
through  the  interaction  terms  now  present  in  the  characteristic  equation  of  the 
overall  system.     This  equation  should  be  investigated  for  zeroes  in  the  right 
half  complex  plane. 

4.  System  errors  due  to  disturbances  entering  the  system  can  be  large. 
However,  these  large  errors  can  be  reduced  considerably  by  addition  of  a 
suitable  high  gain  loop. 

5.  The  system  error  equation  indicates  that  the  overall  system  can  be 
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stable  for  command  inputs,  x.,  and  yet  be  unstable  for  the  disturbance  inputs, 
(X^and  /\J/j  .    In  any  event,  a  complete  error  analysis  for  disturbance  inputs 
should  be  made. 

6„         Zeroes  in  the  right  half  p  plane  of  the  co-factors  of  <£";;  can  create  compen- 
sation problems  through  the  open  loop  relation,  d'Ll  =  z±L  .    Naturally,  c.  must  be 
chosen  so  that  there  are  no  zeroes  in  the  right  half  p  plane  of  1  +■  d...    Making 

the  proper  choice  of  c     can  be  a  formidable  problem  in  some  instances  with 

ii 

the  result  that  no  practical  compensation  can  be  found. 

7.        Interchanging  the  rows  of  the  control  vector,  6    ,  will  alleviate  the 

impractical  compensation  problems  in  some  cases. 

8o         When  compensation  problems  are  encountered  and  changing  the  rows  of 

the    b    vector  does  not  alleviate  the  problem,  the  non-interaction  theory  can 

not  be  used  for  the  design  of  that  particular  "channel"  or  open  loop  part  of  the 

system  which  is  causing  the  problem,,    Some  other  approach  must  be  used  to 

make  the  "channel"  closed  loop  stable.    Normally,  the  end  result  of  deviation 

from  the  non-interaction  theory  is  system  interaction,  which,  it  is  hoped,  can 

be  made  smalL 
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APPENDIX    I 


A.   Linearize   Longitudinal   Equations 

The  actions  for  longitudinal  motion  in  the  general 
or   large  motion  case   ln  ,   ^^   ^^   ^  ^ 

nave   been  shown   in  P.ef.   5"    to   be; 

fr-3)      z  rwy  =  o  =  -i^  q  +  [r„  -r«x]ftp  +.  r,x  D*fc-  f*j  4-w  *r 

Considering  small  perturbations  froL  a  steady  state 
ongitudinal  flight  oondition.  the  above  aquations  can  b, 
linearized  as  follows: 

1.  Assume  the  following  conditions: 

2.  Assume  tne  following  steady  Itate  plus  perturbation 
relations  for  the  system  variables: 

Vp=Vp,^      j   Q„,<^ 
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where  ^P(^{  6^/e  JT 
variables.     u    " 


are  the  perturbation 


3.  The  dimensionless  coefficients,  Xs,Zs,Yg,Ns,  Ls ,  and 
Ms  are  functions  of  o<  v*'  l3  &  (Js  and  R_.  By  using 
non-dimensional  angular  rates, ^.P;  Q   and  R  , 
then  Xs ,  Zs  and  Ma  become  functions  of  &(    6(      Q, 
and  control  inputs  only.  Thus,  neglecting  the  lesser 
important  derivitives,  the  dimensionless  coefficients 
become : 

4.  In  the  wind  axis  system,  it  is  true  that: 
(r-7)    Qw   =  Q^ 

(X-8)    <k     =.  Q-^ 

5.  Linearizing  equations  (3>l)  through  (x~3)  about 
the  steady  state  condition  (subscript  one)  and 
rearranging  in  operator  form  with  D = d/dt  provides: 

(X-9)  bu  ($  4-  b^ff  4-bta  §^  -^b,^  q  =  J„  lr   +  4  4 
*  a)  herd  i 
b„  -    ~  -  kw  - 


^j 


j^-n^-iicos^-e^fGj 


3 
ft 


2no 


3    J 


b,*  =  o 


v&w 
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=   jL  l/^)  cos*',  #)5,^  1 


'w «  ^^  = 


4*  =    k^e  = 


4lt),^-{§J,— : 


e^  c*. 


co-Rju_c: 


bN  =  -  ^  .  -  [^£»«  ^  **%  &.J. 

b*,  =    o 

b34-   =     ^"^     = 


;X1J 


-p-  eyp>Sc  c 


Mq 


£3I    ^  o 


P^Sc  Cm 


o2T, 


« 
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The  fourth  equation  is  found  from  equation  (J--7) 
and  (X-S)  as  follows: 

c\v  —  Q  —  oc  which  when  linearized 

with  6U/  =  o/   -  Q  ="  O  provides: 

Rearranging  and  using  operational  form  results  in: 

*•  —      — — 

(X-  \Z)     b4,  /0j>    4-b4g  Ot    -h  b-43  BuJ  ■f-  fc*M-  £   =    ^4|  <?r+"  ^42  ^€ 

b^<  —  o 

b*z  =  T>    s  ^/df- 

b43  =  D     2  d/o/t- 

h*4  =■  -I 

■f-41      =    O 

B.    Determination  of   the   Elements   of   the   Aircraft  Matrix 

The  vector-matrix  equation  for   controlling  airspeed, 
AJp  ,   and  the  glide   path  angle,  Sw  ,   using   power   and   elevator 
inputs,   6T   and    £e    is: 


(X-,3) 


(I)  - 


*J 


or 


y 


=lA]S 


V/hen  the  relation  above   is   expanded,   the   result    is: 
(XT-I4-)         ^Tp    =    #u  ^t    +■  £•*  <5e 

(r-is)        8^  =  dxx  £t    +  &^  S& 

The  elements,  #lj  ,  in  equations  (T-14)  and  (x"-/5)  are 
found  by  solving  equations  (X-S)  through  (X-izJ  simult- 
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aneously  f  or  /l4  and  Gw  .  Using  determinant  methods  yields 

ft*/  =       k|    b3t      -Ga.51  t>34 

0      W*         3  _| 


/\Jp  c 


"^e  b3i      0       b34 

0  k^i     b4*    -I 


A 


A 


where 


A  = 


Ui 

Jfct 

1 

O 

bll 

biz. 

b23 

b^4 

t>3l 

baa 

0 

b3+ 

o 

b42. 

bi+3 

-1 

Expanding  and  simplifying  the  notation  results  in: 


A  A' 

(X"'7)    &.,,  f>,.te)-UFJj/T  +[U&)-Up)K 

A'  A' 

(A)  -  [b«  b34  b4z  4-  U4  b^  b43  f  b52  bz3  -  b43  b34  t«J 

(B)  =  Lbi4  ^  +"  baa.--.ba4.U3  b,  J 

(C)    =[bz4U2   4- b^^b.^  ba-U3b24b.il 

(E  )    =  [biA  b=n  Ui  -^b^  bji  -  bw  ba-,,  -  b4i  b34  k. 

(  F)   =  [_  Ui  U_  -  b„  b32.  -  b44  b34  b„J 

(&)  =[l>uki-b>,  Ua.-U1.k4  b,,] 
A'  -[K(A)~U,(B)4.b3,(C)j 
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From  equations  (1-16)  and  (1-17)  it  is  clear  that  the 
matrix  elements  are: 

ft    .  J-„(A)  -£,(&) 

(1-18)    "•"         ^ 

(i-i9)     d  =  -t«(£lr_ki§] 

A' 
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APPENDIX   II 

CALCULATION  OF  THE   F-86   STABILITY  DERIVITIVES  AND 
DETERMINATION  OF   THE   COEFFICIENTS,    bjj   AND   t±* ,OF  APPENDIX  I 

A.  Physical   characteristics   of  the   F-86. 

The   following  data    has  been  obtained  from  Ref .  6   . 

S  =   268  ft? 

W  =   15,500   lbs.    (lightly   loaded) 

b  =  37.1   ft. 

c  =  e.09    ft. 

I  =  86E5    slug-ft? 

I  -  29400    slug-ft? 

Izx  ~  2700   slug-ft? 
I2Z  -  35300   slug-ft? 

B.  Assumed  flight  condition. 

Altitude  =  3ea  level 

Mach  No.  =  .246 

density  =  .002378  slug/ft3 

V  -  163  knots  =  275  ft/sec 

Pj.  rx=  f(r  ^  =.  i±  =  *,  =.  e0Jt  ^o 

©oo  =   —  3     =r   —  .OS*L3  radians 

The  above  data  is  for  an  F-86  proceeding  down  a  three 
degree  glide  slope  at  163  knots  in  a  "clean"  configuration. 

C.  Lift  and  Drag  Coefficients  at  Trim;  Stability  Derivitives 
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2.  From  the  drag  polar  using  C/  =  .-577 as  the  entering 

gument:    u)    —  .0445 
vim 


Therefore,     Drag  ss.  1152  lbs. 

erf,  =  9.5° 
3.  Stability  Derivitives 

The  following  stability  derivitives  were  obtained 
directly  or  indirectly  from  Kef .  6  . 

0^=1.065  per  degree   =  3.72  per  radian 
Cn  /  =  fy*    ^   ^-    —   .361  per  radian 
Caa^j   -  <*c*    ^^     =  -.433  per  radian 


Cnq—"^'^   P«r  radian 

^?£^^     ~     -51  Per  radian 

££<:  =5.1   per   radian 
C/n^=-1.08  per  radian 

CxZt  =•    ,55   Por  ra<*ian 

C/_  =—.1775   per   radian 

CMs   h   C"£e  =   -.945  per  radian 

C/   -  C/c  -  -CmS*   ««356  per  radian 

4.  Coe  fficients  of  longitudinal  equations 

\~T~he:   cedd  t-fon*-l    c*.3£njrrtf>+/oin£     Ore   :\ 

I  *,  -  Q.  -  ^  -  o  J 


7? 
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b-  fa*  e&c^-R* 


n^i    _ 


zn 


n 


-  ffi-CoJ-   ^ 


K 


VoS   —     — ,  foO&>~7  per  <tfad.t a ru 


C"     Avr.  .  1 


'V/V 


to 


<2-  Srsv    = 


3)g 


I  I      I 


»1 


*  LifMGfcrB,  **]-  £f£^)_  W, 


—  .709    -L 


7?'\ 


Vn,  ^ 


/H2*> 


-  (V,  )  =    -   0147  - 


nxd. 


j'   ^S 


'Wpt' 


e    — 


>6 


lb.  K^v 


2. 


> 


=  e^iff:  c«.  =  ~.jj3  x 


7n.  A^o/ 


^xyy 


/vj^ 


/U*o 


V.  K< 


n 


z^^' 


^x: 


-/V7 


a  £ 


/r 


* 


=     —  .^67    — 


/4^C 


jo€ 
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D.  Coefficients  of  Equations  (1-9)  through  (1-12). 


bu-l 

D/6   - 

•     Kvv    -     .0311 

b12=" 

"^  ; 

=      .6067 

b13  = 

^Sw 

=  1.00 

b14  = 

0 

•B. 

-    .318 

b21  = 

-Kzv 

=    .00085 

b22* 

D-  K^s  D  4-   .709 

b23~" 

=  -.00613 

b24^' 

L. 

Kzq]=-.9853 

f21  = 

zST 

=  ♦.00624 

f  22  = 

Kzse 

=  -.07 

b31  =  ' 

_Kqv 

=  0.00 

b32  =  - 

"[DKq*  ■*■    KqJ  =  -113 

b33  = 

0 

*>34* 

D-  K 

=     D    4-  .567 

f31- 

0 

f32- 

KQS<l 

=,-6.74 

b41- 

0 

b42  =* 

D    s 

*Mt 

b43   = 

.D    a 

dUt 

b44    = 

0 

f41    = 

0 

f42   = 

.0 
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APPENDIX  III 


ELEMENTS  OF  THE  SERVO  AND  TRANSDUCER  MATRICES,  T  and  S 


A.  Servo  Matrix  Els  merits,  s^. 

1.  The  relation  between  elevator  deflection,  X  ,  and  the 

signal  to  the  elevator  servo,  <3  ,  is  defined  as  s   . 

A  typical  elevator  servo  system  can  be  described  as  a 

second  order  system  with  natural  f requency ,  U)n    ,  and  damping 

-p 
ratio,  jg  ,  The  phase  relation  for  such  a  system  is  given 


(III-l)    0  =-+«•"■' |  |_  u£ 


If  Jq  ~  *5»    ttL0  closed   loop  transfer   function  has   a   db 
versus   log  w   characteristic  a3   follows: 


maani4u<te 


A  co 


— ♦, 

It  follows  that  if  the  elevator  servo  is  designed  such  that 
the  natural  frequency,  tJne  ,  is  larger  than  the  expected 
range  of  system  frequencies  ,  Auj  ,  or  u)  «.  U)nethen  equation 
(III-l)  can  be  approximated  by 

(in-a)    0=-W'[-%^ 

If  the  damping  ratio  is  taken  as  .5?then  the  phase  relation 
becomes : 


(III-3) 


(26  -  -4cld' 


_i   U)   __ 


u)n. 


=  "fan"1  ^Koo 


inhere   O^—  — 
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It  may  be  shown  that  the  phase  relation  of  (III-3)  is 
that  phase  relation  for  a  simple  first  order  system  of  the 
form:  I  / 

(III-4)      P     ^T+~[  T^! 

In  conclusion,  if  the  damping  ratio  is  chosen  as  .5  , 
and  the  natural  frequency  of  the  servo  is  taken  as  25  rad/sec 
thenvthe  expected  range  of  system  frequencies  the  second 
order  elevator  servo  system  may  be  approximated  by  the  first 
order  system:  ^ 

(III.5)  s/p) « 4- =  — J — 

2.  The  relation  between  effective  throttle  perturbation  or 
thrust,  dj-   ,  and  the  actual  throttle  perturbation,  £T   ,  is 
simply  described  by  the  following  first  order  system: 

(III_6,    fl  = 

^t        rp-f-i 

where  'Y-   the  time  delay  in  obtaining  the  new 
engine  EFfcl  which  is  representative  df  the 
thrust. 

The  above  relation  is  pointed  out  in  Ref .  7  and  is 

generally  true  for  most  turbojets.  In  addition,  a  typical 

time  delay  of  two  seconds,  T"—  2.0  ,  was  obtained  from  this 

reference.  Thus: 

(III-7!   -  -  -  -  £  _     i 


Vp)  =  64- 


B.  Transducer  Matrix  Elements,  t**  • 

1.  The  transducer  that  relates  actual  indicated  aircraft 


G4 


velocity,  /\J^  ,  to  the  measured  aircraft  indicated  velocity 
is  given  by: 

(III-6)    ^|/p)=   "=£-  - 


aJ7>  >2y+  I 

Such  a  transfer  function  assumes  no  position  errors  in  the 
airspeed  indicator  nor  any  dynamic  errors  in  the  instrument. 
The  time  "delay,  '?  =  .2  sec  ,  is  due  to  the  pressure  trans- 
mission delay  in  the  pitot  system. 

The  basis  of  the  foregoing  statements  is  derived  from 
a  discussion  of  airspeed  indicators  and  measuring  devices 
contained  in  Ref .  2   • 

2.  The  transducer  that  relates  actual  glide  path  angle, 
B^  ,  to  measured  glide  path  angle,  G^   ,  is  defined  as  t22# 
It  is  assumed  that  the  angle  G     and  the  angle  of  attack,  o( 
are  obtained  without  error  from  a  gyro  and  an  angle  of 
attack  indicator.  These  two  quantities  are  added  together 
to  produce  0^ .  Therefore: 

{iii-9)     if)     a. 

0to 
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